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Electron-positron annihilation into hadrons plus an energetic photon from initial-state radiation allows the 
hadronic cross-section to be measured over a wide range of energies at high luminosity meson factories. Weighted 
integrals over this cross-section are a decisive input for electroweak precision tests. A Monte Carlo event generator 
called PHOKHARA has been developed, which simulates e + e~ — * hadrons + photon(s) at the NLO accuracy. 
The latest tests and upgrades are presented in this paper. 



1. INTRODUCTION 

The total cross-section for electron-positron 
annihilation into hadrons is one of the fundamen- 
tal observables in particle physics. Its high energy 
behaviour provides one of the first and still most 
convincing arguments for the point-like nature of 
quarks. Its normalization was evidence for the 
existence of quarks of three different colours, and 
the recent precise measurements even allow for an 
excellent determination of the strong coupling at 
very high and intermediate energies (e.g. || 
and refs. therein) through the influence of QCD 
corrections. 

Weighted integrals over this cross-section with 
properly chosen kernels are, furthermore, a deci- 
sive input for electroweak precision tests. This 
applies, for example, to the electromagnetic cou- 
pling at higher energies or to the anomalous mag- 
netic moment of the muon. 

Of remarkable importance for these two appli- 
cations is the low energy region, say from thresh- 
old up to centre-of-mass system (cms) energies 
of approximately 3 GeV and 10 GeV, respec- 
tively. Recent measurements based on energy 
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scans between 2 and 5 GeV have improved the 
accuracy in part of this range. However, similar, 
or even further improvements below 2 GeV would 
be highly welcome. The region between 1.4 GeV 
and 2 GeV, in particular, is poorly studied and no 
collider will cover this region in the near future. 
Improvements or even an independent cross-check 
of the precise measurements of the pion form fac- 
tor in the low energy region by the CMD2 and 
DM2 collaborations would be extremely useful, in 
particular in view of the disagreement found with 
the analysis based on isospin-breaking-corrected 
r decays 0j . Since this dominates in the analysis 
of the muon anomalous magnetic moment, this 
kind of studies will help to clarify the situation 
with respect to the experimental measurement 
of this quantity. 



2. RADIATIVE RETURN AT MESON 
FACTORIES 

Experiments at present electron-positron col- 
liders operate mostly at fixed energies, albeit with 
enormous luminosity, with BABAR and BELLE 
at 10.6 GeV, CLEO-C in the region between 3 and 
5 GeV, and KLOE at 1.02 GeV as most promi- 
nent examples. 

This peculiar feature allows the use of the ra- 



2 



diative return, i.e. the reaction 

e + (px)+e - (pa) -» 7(ki)+7*(<9)H hadrons) , (1) 

to explore a wide range of Q 2 in a single experi- 
ment §gg§|@. 

Nominally an invariant mass of the hadronic 
system between 2m.Tr and the cms energy of the 
experiment is accessible. In practice, to clearly 
identify the reaction, it is useful to consider only 
events with a hard photon — tagged or untagged 
— which lowers the energy significantly. 

The study of events with photons emitted un- 
der both large and small angles, and thus at a 
significantly enhanced rate, is particularly attrac- 
tive for the 7r + 7T~ final state with its clear sig- 
nature, an investigation performed at present at 
DA$NE 0,||,^30 ■ Events with a tagged pho- 
ton, emitted under a large angle with respect to 
the beam, have a clear signature and are thus par- 
ticularly suited to the analysis of hadronic final 
states of higher multiplicity p5l . 

3. MONTE CARLO SIMULATION 

To arrive at reliable predictions including kinc- 
matical cuts as employed by realistic experiments, 
a Monte Carlo generator is indispensable. The 
inclusion of radiative corrections in the generator 
and the analysis is essential for the precise ex- 
traction of the cross-section. For hadronic states 
with invariant masses below 2 or even 3 GeV, it 
is desirable to simulate the individual exclusive 
channels with two, three, and up to six mesons, 
i.e. pions, kaons, etas, etc., which requires a fairly 
detailed parametrization of the various form fac- 
tors. 

A first program, called EVA, was constructed 
some time ago Q to simulate the production 
of a pair of pions together with a hard pho- 
ton. It includes initial-state radiation (ISR), 
final-state radiation (FSR), their interference, 
and the dominant radiative corrections from addi- 
tional collinear radiation through structure func- 
tion (SF) techniques [fl6| . This project was con- 
tinued with the construction of a generator for 
the radiative production of four pions . 

As a further development of this project a new 
Monte Carlo generator called PHOKHARA @ 




Figure 1. One-loop corrections to initial-state ra- 
diation in e + e~ — > 7+ hadrons. 



has been constructed, which includes, in contrast 
to the former generators, the complete next-to- 
leading order (NLO) radiative corrections. The 
first version of PHOKHARA incorporates ISR 
only and is limited to -k + tt~^{^) and fi + fi~"f('j) 
as final states. PHOKHARA exhibits a modu- 
lar structure that simplifies the implementation 
of additional hadronic modes or the replacement 
of the currents(s) of the existing modes. 

In this paper we show how PHOKHARA is con- 
structed and present some results obtained with 
its current version, as well as a comparison with 
the other aforementioned generators. Finally, the 
new features that will be part of the next ver- 
sion of PHOKHARA are outlined and some pre- 
liminary results are presented. Further details 
will be given, however, in a forthcoming publi- 
cation |l8| . These programs and the future ver- 
sions of PHOKHARA can be downloaded from 
tittp : //cern. ch/german. rodrigo/phokhara/. 



4. NLO CORRECTIONS TO ISR 

At NLO, the e + e~ annihilation process ([!]), 
where the virtual photon converts into a hadronic 
final state, J*(Q) — > hadrons, and the real one is 
emitted from the initial state, receives contribu- 
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Figure 2. Emission of two real photons from the 
initial state in e + e~ annihilation into hadrons. 



tions from one-loop corrections (see Fig. [l]) and 
from the emission of a second real photon (see 
Fig. 1). 

After renormalization the one-loop matrix el- 
ements still contain infrared divergences. These 
are cancelled by adding the contribution where 
a second photon has been emitted from the ini- 
tial state. This rate is integrated analytically in 
phase space up to an energy cutoff < w^/s 
far below t/s. The sum is finite; however, it de- 
pends on this soft photon cutoff. The contribu- 
tion from the emission of a second photon with 
energy E 7 > w^fs completes the calculation and 
cancels this dependence. 

In order to facilitate the extension of the Monte 
Carlo simulation to different hadronic exclusive 
channels, the differential rate for the virtual and 
soft QED corrections is cast into the product of 
a leptonic and a hadronic tensor and the corre- 
sponding factorizcd phase space: 



da 



2s 



L^ v d$ 2 (pi,P2;Q,ki) 



H" v d$„(Q;<?i,-,?») 



dQ 2 
2n ' 



(2) 



where d<fr n (Q\qi,-,q n ) denotes the hadronic n- 
body phase space including all statistical factors 
and Q 2 is the invariant mass of the hadronic sys- 
tem. 

The physics of the hadronic system, whose de- 
scription is model-dependent, enters only through 



the hadronic tensor 

= JfJ u+ , (3) 

where the hadronic current has to be 
parametrized through form factors |?| JTg| , pc| , pl| . 

The leptonic tensor, which describes the NLO 
virtual and soft QED corrections to initial-state 
radiation in e + e~ annihilation, has the following 
general form: 

(47ra) 2 
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Terms proportional to are absent as a conse- 
quence of current conservation. The scalar coeffi- 
cients aij and a_i allow the following expansion 



(0) 



where af^ provide the leading order (LO) con- 
tribution and the imaginary antisymmetric piece 
proportional to a_i appears for the first time at 
NLO. 

As an alternative one can replace the Cartesian 
basis (||) by a basis derived from the three circular 
polarization vectors of the virtual photon e l and 
e±: 

(Ana) 2 v 



a (i) 
— a) 

7T 



I J 



a (i) 
a-x = — a_{ , 

7T 



(5) 



Hj £ i E j 



i,j = L,± , (6) 



where only four of the scalar coefficients are in- 
dependent 



a-L- = a L + 



a-L = a+L 



a = a++ , a |_ = a+_ . 

An expansion similar to (||) holds for these scalar 
coefficients. The relationship between the com- 
ponents in both bases as well as expressions for 
these coefficients can be found in f2^Jj3|| . The 
trace of the leptonic tensor, which is related to 
the cross-section after angular averaging of the 
hadronic tensor, is particularly simple in the sec- 
ond case: 



(Ana) 2 



(a LL + 2a ++ ) .(7) 
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Note that the imaginary part of L^, which is 
present in the coefficients or a_i only, is of in- 
terest for those cases where the hadronic current 
receives contributions from different amplitudes 
with non-trivial relative phases. This is possible, 
e.g. for final states with three or more mesons or 
for pp production. 

The matrix elements for the emission from the 
initial state of two real hard photons, i.e. E 1 > 

e + (pi) + e-(p 2 )^7*(Q)+7(fci)+7(fe) , (8) 

are calculated numerically following the helicity- 
amplitude method with the conventions intro- 
duced in HJ.^H. As a test, the square of this 
matrix element averaged over initial particle po- 
larization has also been calculated using the 
standard trace technique and tested numerically 
against the helicity method result. 

The virtual plus soft contribution and the hard 
one depend separately on the soft photon cut- 
off w used to regulate the infrared divergences of 
the virtual diagrams. The former shows a log- 
arithmic w dependence. The second, after nu- 
merical integration in phase space, exhibits the 
same behaviour, whereas their sum must be in- 
dependent of w. To explicitly demonstrate this 
w-independence is therefore a basic test of the 
program. Then the value of w that optimizes the 
event generation, avoiding at the same time the 
appearance of unphysical negative weights, is de- 
termined. 

Table |l| presents the total cross-section for ra- 
diative production of a pair of pions calculated for 
several values of the soft photon cutoff at three 
different cms energies for the kinematical cuts 
from Table ^. The excellent agreement, within 
the error of the numerical integration, confirms 
the uj-indcpcndcncc of the result. A value around 
w = f0~ 4 seems to be the best choice [flTll . 
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5. LL VERSUS NLO 

The original and default version of EVA §], 
simulating the process e + e~ — > 7r + 7r"7 at LO, al- 
lowed for additional initial-state radiation of soft 
and collinear photons by the structure-function 
method p6l. By convoluting the Born cross- 



Figure 3. Comparison between the collinear ap- 
proximation by structure functions and the fixed- 
order NLO result for different cms energies. Cuts 
from Table @. 
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Table 1 

Total cross-section (nb) for the process e + e~ — > 7r + 7r~7 at NLO for different values of the soft photon 
cutoff. Only initial-state radiation. Cuts from Table |[ 

w =L02 GeV 4 GeV 10.6 GeV 

IST 3 2.0324 (4) 0.12524 (5) 0.010564 (4) 

10~ 4 2.0332 (5) 0.12526 (5) 0.010565 (4) 

10~ 5 2.0333 (5) 0.12527 (5) 0.010565 (5) 



Table 2 

Kinematical cuts applied at different cms energies: minimal energy of the tagged photon (£™ m ), angular 
cuts on the tagged photon (# 7 ) and the pions (9^), and minimal invariant mass of the hadrons plus the 
tagged photon (M%+ ) 
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Table 3 

Total cross-section (nb) for the process e + e~ — * ir + ir~j at LO, NLO and in the collinear approximation 
via structure functions (SF), with the cuts from Table Q. Only initial-state radiation. NLO(2) gives the 
NLO result with the same cuts as NLO(l) but for the hadron-photon invariant mass unbounded. 





y/s =1.02 GeV 


4 GeV 


10.6 GeV 


Born 


2.1361 (4) 


0.12979 (3) 


0.011350 (3) 


SF 


2.0192 (4) 


0.12439 (5) 


0.010526 (3) 


NLO (1) 


2.0332 (5) 


0.12526 (5) 


0.010565 (4) 


NLO (2) 


2.4126 (7) 


0.14891 (9) 


0.012158 (9) 



Table 4 

Total cross-section (nb) for initial-state radiation in the process e + e~ — > /i + /i _ 7 at LO, NLO (1) and 

NLO (2) with the cuts from Table |, the pions being replaced by muons. 

y/s =1.02 GeV 4 GeV 10.6 GeV 

Born 0.8243(5) 0.4690(6) 0.003088(6) 

NLO (1) 0.7587(5) 0.4449(6) 0.002865(6) 
NLO (2) 0.8338(7) 0.4874(14) 0.00321(6) 
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Table 5 

Total cross-section (nb) for the process e + e~ — * 
7T+7T-7 at yfs =1.02 GeV in NLO and in the 
collinear approximation (SF) as a function of 
the cut on the invariant mass of the hadron + 
tagged photon M 2 +w _^. Only initial-state radia- 
tion. Minimal energy of the tagged photon and 
angular cuts from Table ||. 
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section with a given SF distribution, soft pho- 
tons are resummed to all orders in perturbation 
theory and large logarithms of collinear origin, 
L = log(s/?7ig), are taken into account up to 
two-loop approximation. The NLO result, be- 
ing a complete one- loop result, contains these 
logarithms in order a and additional subleading 
terms, which of course are not taken into account 
within the SF method. Although small, these 
subleading contributions become important for 
measurements aiming at an accuracy better than 
1%. Furthermore, a NLO generator is more suit- 
able for comparison with an experimental setup. 
While in the SF approach the extra collinear pho- 
ton emission is integrated out and some informa- 
tion is lost, in a fixed-order calculation the full an- 
gular dependence is kept and energy-momentum 
is conserved, i.e. the sum of the momenta of the 
generated outgoing particles agrees with the in- 
coming e + e~ momentum. 

Table || presents the total cross-section for 
e + e~ — > 7r + 7r~7(7) calculated at LO and NLO for 
three different cms energies with the kinematical 
cuts listed in Table ||. Two NLO predictions are 



shown. The first one, NLO(l), which can be com- 
pared with the SF result derived from EVA, in- 
cludes a cut on the invariant mass of the hadrons 
plus the tagged photon. The last was introduced 
in |^| in order to reduce the kinematic distortion 
of the events due to the extra collinear emission. 
The second prediction, NLO (2), is obtained with- 
out this cut. The Q 2 dependence of the difference 
between the NLO(l) prediction and the SF result 
is showed in Fig. [| The size and sign of the NLO 
corrections do depend on the particular choice of 
the experimental cuts. Hence only using a Monte 
Carlo event generator can one realistically com- 
pare theoretical predictions with experiment and 
extract R(s) from the data. 

The results of EVA and those denoted NLO(l) 
for the total cross-section are in reasonably good 
agreement. Both of them are clearly sensitive to 
the cut on M 2 + . This cut dependence is dis- 
played in Table ||. Remarkably, the typical differ- 
ence between the results of the two programs is 
clearly less than 0.5% for most of the entries. 

The systematic uncertainty of the program, 
due to inadequate treatment of ISR, namely miss- 
ing leading logarithmic higher order corrections 
and lepton pair production, is conservatively es- 
timated to be of around 0.5% in the total cross- 
section jl7jj . 

6. MUON PAIR PRODUCTION 

Results for the total cross-section of muon pair 
production are listed in Table f|. The radiative 
muon cross-section can be used for a calibration 
of the pion yield. A number of radiative correc- 
tions are expected to cancel in the ratio. For 
this reason we consider the ratio between the 
pion and the muon yields, after dividing the for- 
mer by |F 7r (Q 2 )| 2 (l-4m 2 r /Q 2 ) 3 / 2 , the latter by 

4(1 + 2m 2 /Q 2 )y / l - 4m 2 /Q 2 . In Fig. |a, we con- 
sider the full angular range for pions and muons, 
with 6> 7 between 5° and 21°. Clearly all radia- 
tive corrections and kinematic effects disappear, 
up to statistical fluctuations, in the leading or- 
der program as well as after inclusion of the NLO 
corrections. 

In Fig. Elb an additional cut on pion and muon 
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Figure 4. Ratio between pion and muon yields, after dividing through their respective R-ratio. a: no 
cuts on pion and muon angles; b: with angular cuts on pion and muon angles. 



angles has been imposed. As demonstrated in 
Fig. ^b, the ratio differs from unity once (iden- 
tical) angular cuts are imposed on pions and 
muons, a consequence of their different angular 
distribution. To derive the pion form factor from 
the ratio between pion and muon yields, this ef- 
fect has to be incorporated. However, the cor- 
rection function shown in Fig. ||b is independent 
from the form factor, and hence universal and 
model-independent (ignoring FSR for the mo- 
ment). 

7. UNTAGGED PHOTONS 

Both EVA and PHOKHARA were initially de- 
signed to simulate reactions with tagged photons, 
i.e. at least one photon was required to be emit- 
ted under large angles. The extension of these re- 
sults to untagged photon events, i.e. for photons 
emitted at arbitrary small angles, was recently 
investigated 

An important ingredient in the extension of 
the NLO Monte Carlo program PHOKHARA 
to small photon angles is the evaluation of the 
virtual corrections to reaction ([!]) in the limit 
TOg/s <C 1, which are equally valid for large 
and small angles. Compact results for the one- 
loop two-, three- and four-point functions that 
enter this calculation can be found in the litera- 
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Figure 5. Comparison of the virtual+soft and 
hard contributions to the 7r + 7r~ differential cross- 
section with inclusive analytical results. Soft pho- 
ton cutoff: w = 10 -4 . 
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ture fl27| , p8|l for arbitrary values of m 2 e / s. How- 
ever, the combination of these analytical expres- 
sions with the relevant coefficients is numerically 
unstable in the limit of small mass and angles. 
A compact, numerically stable result, valid for 
an arbitrarily small photon angle, is therefore re- 
quired. As a consequence of the highly singular 
kinematic coefficients, terms proportional to m 2 
and even ra\ must be kept in the expansion, which 
will contribute, after angular integration, to the 
total cross-section even in the limit rn^/s — ► 0. 

Inclusive NLO calculations can be used to test 
the performance of the simulation at small an- 
gles. Figure || shows the comparison of the differ- 
ential cross-section for the tt + tt~ mode with the 
analytical results from 29 3(J. The agreement is 
excellent. 

8. ISR DOMINANCE AND FSR 

Final-state radiation should be regarded as the 
background of the measurement and would be re- 
quired for the complete simulation of the reac- 
tion. The complete leading order matrix element 
squared is given by 

|-M| 2 = |Mi S r| 2 +|X fsr | 2 +2Rc[X isr a4 sr ] .(9) 

FSR and its interference with ISR were already 
included in EVA g for the two-pion case. The 
pions were assumed to be point-like, and scalar 



QED was applied to simulate photon emission off 
the charged pions. There it was demonstrated 
that suitably chosen configurations, namely those 
with hard photons at small angles relative to the 
beam and well separated from the pions, are dom- 
inated by ISR. FSR can therefore be reduced to 
a reasonable limit or be controlled by the simula- 
tion. 

At B-factories, where one has to deal with very 
hard tagged photons, the situation is even bet- 
ter because the kinematic separation between the 
photon and the hadrons becomes very clear. For 
events where hadrons and photon are produced 
mainly back to back the suppression of FSR is 
naturally accomplished. 

The third term in the right-hand side of Eq. (g), 
the ISR-FSR interference, is odd under charge 
conjugation and its contribution vanishes under 
angular integration. It gives rise, however, to 
a relatively large charge asymmetry and, corre- 
spondingly, to a forward-backward asymmetry 



A{9) 



N* + {6) - N* + (ir - 9) 



(10) 



N* + (9) + N* + (tt - 9) ' 

The asymmetry can be used for calibration of the 
FSR amplitude, and fits to the angular distribu- 
tion A{9) can test details of its model dependence. 

This is illustrated in Figs. ^ and [?], where the 
angular distributions of 7r + and fi + respectively 
are shown for different kincmatical cuts. The 
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Figure 8. Diagrams contributing to the 4-7T hadronic current. 



angles are denned with respect to the incoming 
positron. If no angular cut is applied the angu- 
lar distribution in both cases is highly asymmet- 
ric as a consequence of the ISR-FSR interference 
contribution. If cuts suitable to suppress FSR, 
and therefore the ISR-FSR interference, are ap- 
plied the distributions become symmetric. These 
plots have been obtained with the new version 
of PHOKHARA, which now incorporates at LO 
FSR, and its interference with ISR, for two pions 
(point-like) and muons. 

9. NEW HADRONIC CHANNELS 

Because of the modular structure of 
PHOKHARA additional hadronic modes can 
be easily implemented. The four-pion channels 
(2ir + 2w~ and 27r 7r + 7r~), which give the domi- 
nant contribution to the hadronic cross-section 
in the region from 1 to 2 GeV, are a new feature 
of our event generator. 

Isospin invariance relates the amplitudes of 
the e + e~ — * 2tt + 2tt^ and e + e~ — > 2tt°tt + tt^ 
processes and those for r decays into 7r~37r° 
and 7T + 27r~7r°. The description of the four-pion 
hadronic current follows 0,0. The basic build- 
ing blocks of this current are schematically de- 
picted in Fig. H and described in detail in 0. 

Results obtained with PHOKHARA for this 
channel have been compared with the Monte 
Carlo, which simulates the same process at 
LO and includes additional collinear radiation 
through SF techniques. Typically, differences of 
order 1% are found (see Fig. ||), which are of the 
expected size and of the same order as for the 
two-pion final state p7| . Details of the compari- 
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son and of the performed tests will be presented 
in a forthcoming publication Jl8| . 

The implementation of the three-pion mode is 
under way. Other hadronic channels will be like- 
wise considered in the future. 

10. CONCLUSIONS 

The PHOKHARA Monte Carlo event genera- 
tor has been upgraded. It now includes, besides 
the 7r + 7r~ and channels, also 27r + 27r _ and 

27r°7r + 7r - as final states. FSR, and its interfer- 
ence with ISR, have also been included for the 
first two channels. Furthermore, the simulation 
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of events where the photon(s) is emitted under 
small angles with respect to the beam axis has 
been tested numerically. Further upgrades are in 
progress. 
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